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The sequence for adapting a conjugate mathematical model describing external heat exchange in the work 
space and heat exchange and hydrodynamics of the melt in a glass-melting furnace melting tank with horse- 
shoe-shaped flame direction was examined. It was shown that it is sufficient to use the experimental 
dependences of the air heating temperature and temperature of the glass melt at the outlet of the melting tank 
on the furnace efficiency for adapting the model when the basic modeling boundary conditions have been cor¬ 
rectly assigned. Equations were obtained that establish the correlation of the specific heat rate and maximum 
melting temperature and the specific output of glass melt, and they can be used in constructing an algorithm 
for controlling the thermal work of a glass-melting furnace of the design studied. 


In the current conditions of development of the glass in¬ 
dustry, mathematical modeling of glass-melting furnaces 
should be considered an integral part of design methodology 
[1]. Mathematical models that are sufficiently complete and 
adequate can be used to not only evaluate the reliability of 
the basic characteristics of furnaces but also to optimize 
them and to predict their behavior with consideration of the 
probable nature of many operating parameters of real ob¬ 
jects. Predicting the output characteristics of a designed ob¬ 
ject (output, fuel consumption, melting temperature, etc.) in 
the design stage is due to the necessity of ensuring the re¬ 
quirements for the quality of the glass melt, thermal effi¬ 
ciency of the furnace, and operating time between servicings. 
It makes it possible to evaluate the correspondence of the 
furnace design to the performance loads and purposefully se¬ 
lect the optimum operating parameters. 

In the general case, the detail and accuracy of the mathe¬ 
matical model should correspond to the given quality of the 
object. Superfluous details in the model causes an unjustified 
increase in computer time costs and insufficient detail leads 
to an incorrect result. On the whole, the basic requirements 
for mathematical models can be reduced to goodness of fit, 
accuracy, and economy. Economy (computer efficiency) is 
characterized by the costs for computer time and memory of 
the computer system. Goodness of fit is the reflection of the 
selected properties of the object by the model with a given 
accuracy. Accuracy means the degree of correspondence of 
the evaluations of similar properties of the object and model. 
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The goodness of fit of the model to the sample is always 
limited and is a function of the purpose of modeling. Any 
model does not take into account some properties of the ori¬ 
ginal and for this reason is an abstraction of it. The goodness 
of fit of the model is established by testing the basic laws of 
the object region (for example, laws of conservation) and 
comparing the results of modeling of particular variants with 
the known solutions for these variants. Within the framework 
of the principles of reasonable sufficiency and balance of ac¬ 
curacy, the evaluation of the goodness of fit of a mathemati¬ 
cal model in the first stage consists of establishing the physi- 
cality of the results of the calculations [2, 3]. Here physicali- 
ty means the qualitative correspondence of the mechanisms 
revealed in modeling developed in the object region to the 
representations or available experimental data. On the quan¬ 
titative level, the degree of agreement of the calculated quan¬ 
tities with the real parameters of the full-scale object is deter¬ 
mined by the modeling goals. 

We can distinguish two limiting cases of the use of ma¬ 
thematical modeling of glass-melting furnaces. In the first 
one, the effect of some design (process) parameters of the ag¬ 
gregate on its operation must be established. It is evident that 
modeling in this case will boil down to executing a series of 
comparative calculations, and high accuracy of correspon¬ 
dence of the calculated data with the real values is not re¬ 
quired here [4-6]. It is only important for the effect of the 
investigated parameters on the evolution of heat exchange 
and hydrodynamic processes to be reflected in the model. In 
other words, establishing the physicality is the task of evalu- 
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ating the goodness of fit of a mathematical model in the 
given case, and for this reason, this evaluation can frequently 
only be conducted without involving experimental data. 

In the second limiting case, the mathematical model is 
either used as a constituent part of the system for automatic 
control of the object or for revealing “bottlenecks” in this ob¬ 
ject; for example, for determining fuel consumption, the 
maximum temperate of the protective lining, glass-melt tem¬ 
perature at the furnace outlet, etc. Since the working charac¬ 
teristics of the design elements are taken into consideration 
in the calculation, in the given case, higher requirements are 
imposed on the evaluation of the goodness of fit of the math¬ 
ematical model and they can only be satisfied if special ex¬ 
periments are conducted. Adaptation of a mathematical 
model not only implies refining the initial modeling condi¬ 
tions, but also quantitatively evaluating the goodness of fit of 
the basic integral indexes of operation of a furnace of analo¬ 
gous design. Incorporating a furnace efficiency variable is 
the basic condition for executing this procedure. 

As an illustration of the above, consider the sequence of 
adaptation of a conjugate mathematical model of hydrody¬ 
namics and heat exchange in a glass-melting furnace with 
horseshoe-shaped flame direction which is completely cha¬ 
racterized in [2, 3]. We will use the operating indexes of 
glass-melting furnace No. 1 at Orekhovo-Zuevo Glass Com¬ 
pany as the experimental data [1], The results of mathemati¬ 
cal modeling were used in designing this furnace, and the 
difference between it and the model only consisted of the 
area of the melting tank (F t = 102.4 m 2 ) and the lower heat 
value of natural gas (Q{ = 33.446 MJ/m 3 ). In modeling 
these values were set equal to 115.77 m 2 and 35.042 MJ/m 3 
in [2], The experimental data were obtained for the furnace 
output range of 115.0 - 241.1 tons/day. 

Consider the furnace operating parameters used for eva¬ 
luating the goodness of fit of the modeling results. The con¬ 
trol-measurement complex of a modem glass-melting fur¬ 
nace allows determining its output in fused glass melt P 
(tons/day), fuel consumption B (m 3 /day), temperature of the 
inside of the roof, including the maximum value t im , air 
heating temperature in the regenerator t , and weight-aver¬ 
age temperature of the glass melt at the furnace outlet t g . 
These parameters can be used to calculate the specific fuel 
consumption J? sp or the specific heat flow q sp , as well as the 
specific output of glass melt from 1 m 2 of melting tank sur¬ 
face area P sp (ton/(m 2 • day)). It is thus possible to determine 
the dependences that are integral characteristics of the ther¬ 
mal efficiency of a glass-melting furnace: 

S sp = /( p )> m3 /ton; (1) 

9sp = /(^sp )> kJ/kg. (2) 

Equation (1) can be directly used in the algorithm for 
controlling the thermal work of the furnace. Equation (2) is 
the most universal characteristic of operation of a furnace 
generally accepted in world practice. Using it allows not only 


comparing the operating indexes of different furnaces of the 
same design type but also the results of modeling with expe¬ 
rimental data obtained on samples with values of F t and Q{ 
that differ from the model. On the whole, functions (1) and 
(2) can be used for evaluating the goodness of fit of the 
mathematical model to the industrial sample. 

The convergence of the numerical values of Eqs. (1) and 
to a greater degree (2) for the model and sample is an impor¬ 
tant but not the only condition for testing the goodness of fit 
of the results of the calculation. We know that in practical 
glass making, controlling the thermal work of the furnaces is 
based on some “arbitrary” indexes on which the furnace out¬ 
put and glass quality are to some degree dependent. They in¬ 
clude not only the fuel (or heat) consumption but also the 
temperature distribution in the working space of the furnace 
and its (temperature) level [7], It is customarily assumed that 
the conditions required for the successful occurrence of the 
manufacturing process can be created by correctly assigning 
the level and position of the temperature maximum. 

The temperature in the furnace is one of the most com¬ 
plex control parameters. This is first due to the complexity of 
measuring the real temperature and limiting the number of 
control points, which does not give a sufficiently complete 
representation of the temperature field in the furnace. We 
know that the gas medium in the working space of high-tem- 
perature furnaces is characterized by a nonuniform tempera¬ 
ture field. For this reason, the term “furnace temperature” is 
arbitrary in character, and the notion of “effective tempera¬ 
ture” is used for characterizing the thermal state of the work¬ 
ing space. Second, this is related to the low inertness of the 
object regulated and the mutual effect of the temperatures in 
the furnace and the controlling action — fuel or power con¬ 
sumption. 

To increase the stability of operation of the control sys¬ 
tem, slight “coarsening” of the primary signal is useful. This 
implies selecting a more inert control object. The tempera¬ 
ture of the inner surface of the working space lining, prefera¬ 
bly the furnace roof, can serve as such an object. Assigning 
the temperature curve and especially its maximum based on 
the temperature of the inside of the roof produces several 
positive results. First, more accurate measurement of this 
temperature by mounting the case of the thermocouple flush 
with the inner surface of the lining. Second, objective control 
of the maximum roof temperature is ensured, which prevents 
accidents related to its overheating. Third, the roof tempera¬ 
ture is more inert to brief deviations in the furnace operating 
regidations. And finally, and most important, in certain con¬ 
ditions, the maximum roof temperature is satisfactorily cor¬ 
related with other furnace thermotechnical and performance 
parameters. The maximum value of the temperature of the 
inner surface of the roof and its position over the length of 
the furnace can thus be included among the basic parameters 
in testing the goodness of fit of the mathematical model to 
the industrial sample [7]. 

Elaboration of a mathematical model of a glass-melting 
furnace that satisfactorily fits the sample in a wide output 
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range (for example, 100 - 300 tons/day) presupposes correct 
definition of some of the most important initial conditions of 
modeling: input of heat with the heated air and loss of heat 
with the working glass-melt stream. At constant output 
(P = const), these items can be determined from calculating 
the heat balance of the melting part of the furnace. For vari¬ 
able furnace output, formalization of the dependences 
t a =f(P) and t =f(P) implies adaptation of the model us¬ 
ing the experimental data. The necessity of executing this 
procedure can be seen with the data in Fig. 1, where the cal¬ 
culated values for different conditions of formalizing 
t a =f(P) and t gm =f(P) are shown together with the experi¬ 
mental dependence of 5™. We note that the experimental 

data are approximated by the equation (standard error of ap¬ 
proximation R 2 = 0.9778): 

S s e p x = 239 - 0.6346P + 0.0009P 2 , m 3 /ton. 

The heat expended for glass formation Q[ and heating 
the glass melt in the tank Q" to 1500°C at different furnace 
output was calculated with the equations: 

Q[ =-0.240668 + 11.280498P, kJ; 

Q[= 1.393369 + 27.49743 IP, kJ. 

The method of distribution of Q[ and Q" in the tank is 
given in [2], 

If the mathematical model is restricted by the first limit¬ 
ing case of its application, then the simplest method of defin¬ 
ing the initial conditions with t a =/(P) and t„ m =f(P) is 
possible. It consists of the fact that these parameters are de¬ 
fined in the form of fixed quantities whose concrete values 
should be correlated (from an experiment) with the furnace 
output used. If P = var and varies in a relatively wide range, 
then implementation of this variant of defining the initial 
conditions, for example, t a = 1200°C and t gm = 1300°C (see 

Fig. 1, curve 2) will lead to an important difference in the 
calculated and experimental data, both qualitatively and 
quantitatively. 

The source of the calculation error (maximum deviation 
of 9.45%) is primarily due to the fact that at low furnace out¬ 
put, the amount of heat introduced in the regenerator by the 
products of combustion cannot ensure the given level of 
heating of the air in combustion. In addition, it is also prede¬ 
termined by the iteration method of determining the fuel con¬ 
sumption at which the calculations terminated when the 
weight-average temperature of the glass melt in the neck at¬ 
tained the assigned value. 

Another possible method of defining the air heating tem¬ 
perature that to some degree correlates t with furnace output 
and fuel consumption consists of using a relation of the type 

r a = r d -At. 

The temperature of the outgoing fuel combustion pro¬ 
ducts t d is determined by calculation. Quantity At (experi- 


B , m 3 /ton 



P , tons/day 

Fig. 1 . Effect of furnace output on fuel rate: 1 and 2-4 ) experi¬ 
mental and calculated data. 


mental) can be varied within the limits of 150 - 300°C. The 
variant of the calculation for A t= 300°C and t gm = 1300°C 

is shown in Fig. 1 (see curve 3 ). The results of the calcula¬ 
tion indicate that the character of the dependence B sp =/(P) 
did not undergo any major changes in comparison to the pre¬ 
ceding variant (see curve 2 ). 

On the whole, we note that the examined methods of as¬ 
signing the initial conditions with t a and t g m can be used for 
solving a number of problems related to the first limiting 
case of use of mathematical modeling. The value of the cal¬ 
culation error and character of the change in P sp =/(P) suffi¬ 
ciently correctly reflects the real operating indexes of 
glass-melting furnaces. At the same time, going to problems 
of the second limiting case requires more complete both 
quantitative and qualitative correspondence of the calculated 
and experimental data. Attaining this goal is seen in adapta¬ 
tion of the model using the experimental dependences 
t a = f(P ) and t g m =f(P ). 

The analysis of the results of operating an industrial sam¬ 
ple showed that the temperature of heating the air for com¬ 
bustion and the weight-average temperature of the glass melt 
at the furnace outlet are in certain dependences on its output 
(Fig. 2). These dependences can be sufficiently correctly ap¬ 
proximated by second-degree polynomials. The equations 
are as follows for t a (R 2 = 0.9833) and t g m (R 2 = 0.9937): 

t a = 1009.5 + 124.4P- 3.8933P 2 , °C; (3) 

t g m = 1414.6 + 0.3299P-0.0003P 2 , °C. (4) 

The use of Eqs. (3) and (4) as the initial modeling condi¬ 
tions allows obtaining a new calculation expression for func¬ 
tion B sp =f(P) (see Fig. 1, curve 4): 

B sp = 108.223 + 6317617 , m 3 /ton. (5) 
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t, °c 



P, tons/day 

Fig. 2. Experimental dependence of the air heating temperature ( 1 ) 
and glass melt temperature at the neck outlet (2) on furnace output. 


The character of function (5), adapted with experimental 
data, more accurately reflects the change in the fuel value as 
a function of the furnace output (see Fig. 1, curve 1 ). The 
maximum calculation error is 8.68% and the average devia¬ 
tion is equal to 6%. 

We note that curves 1 and 4 (see Fig. 1) were obtained 
for different values of the melting tank surface area and low¬ 
est operating efficiency of the fuel. For this reason, a correct 
comparison of the calculated and experimental data implies 
their representation in the form of dependences described by 
equations of the type of (2). The data in Fig. 3a indicate that 
in generalized form, the results of calculating q =f(P sp ) is 
satisfactorily correlated with the experimental data. The 
higher the specific output of the furnace, the better the con¬ 
vergence of the results of the calculation and the experimen¬ 
tal data. This is because the initial conditions on heat losses 
through furnace structural elements were assigned for fur¬ 
nace output of 300 tons/day [2]. The average deviation of 
the calculation is 3.86%, which indicates the good fit of the 
mathematical model to the real process. The calculated curve 
of function q sp =f(P sp ) can be described by the following 
equation (R 2 = 0.9911): 

q sp = 6668.3 - 1316.3P sp + 181.98P 2 , kJ/kg. (6) 

Adapting the mathematical model with the initial condi¬ 
tions of assignment of t a and t„ m thus allows obtaining the 
equation correlating the heat value with the specific furnace 
output, which indicates the suitability of the model for solv¬ 
ing problems related to the second limiting case. Equation (6) 
can be used in the algorithm of the equation for the thermal 
work of furnaces of the given design in a wide output range. 

The dependence of the maximum roof temperature on 
the specific furnace output is shown in Fig. 3 b. The calcu¬ 
lated values were obtained with consideration of adaptation 
of the mathematical model using the experimental data in 


? sp , kJ/kg 



t, °C 



P sp , tons/(m 2 ■ day) 

Fig. 3. Heat value ( a ) of the maximum roof temperature (b ) on unit 
furnace output: 1 and 2) experimental and calculated data. 


Fig. 2 and approximated equations (3) and (4). Processing 
the results of the calculation showed that function t rm =f{P sp ) 
can be described by the expression (R 2 = 0.9978): 

t lm = 1443.4 - 74.557P sp + 5.67593P S 2 . (7) 

The data in Fig. 3 b indicate that the results of the calcula¬ 
tion correlate satisfactorily with the experimental values both 
qualitatively and quantitatively. The average error of calcu¬ 
lating the maximum roof temperature does not exceed 1.4%. 

The character of the change in function t im =f(P s p ) cor¬ 
responds to the developed concepts not only concerning the 
level of the maximum roof temperature but also the change 
in it as a function of the specific output of modem glass¬ 
melting furnaces with horseshoe-shaped flame direction. For 
example, for P sp = 1.0, 1.5, 2.0, 2.5, and 3.0 tons/(m 2 • day), 
with Eq. (7) we obtain / rm = 1512, 1542, 1569, 1594, and 
1615°C. Two conclusions derive from the reported data 
which indicate the goodness of fit of Eq. (7) to glass-melting 
practice. The first is that the increase in the specific furnace 
output will outstrip (in percentage) the required increase in 
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the melting temperature. The second is that for a flame fur¬ 
nace (without using additional electric heating), the specific 
glass-melt output of 2.5 tons/(m 2 • day) will be the limiting 
value with respect to the temperature of its upper structure. 

In conclusion, it should be noted that the conjugated 
mathematical model of a glass-melting furnace with horse- 
shoe-shaped flame direction developed [2, 3] satisfactorily 
reflects the thermophysical processes that characterize opera¬ 
tion of an industrial sample. In correct assignment of the ba¬ 
sic boundary conditions of modeling (set of initial and 
boundary conditions), it is sufficient to use the experimental 
dependences of the air heating temperature and glass-melt 
temperature at the melting tank outlet on the furnace output 
for adapting the model. 

Equations (5) - (7) can be used in constructing an algo¬ 
rithm for control of the thermal work of a glass-melting fur¬ 
nace of the investigated design. 

The development of systems for controlling the thermal 
work of glass-melting furnaces is correlated with the creation 
of multilevel algorithms. Despite the attractiveness of the 
roof temperature as object of regulation, this parameter does 
not fully reflect all of the complexity of the processes that 
take place in the furnace, and is also not a direct characteris¬ 
tic of any aspects of the quality of the glass melt. Develop¬ 
ment of numerical modeling methods will allow making the 
control algorithm more complex, more reliable, and better fit 
the furnace operating conditions. The possibility of a com¬ 
bined solution to external and internal heat-exchange prob¬ 


lems will form the basis for passing to control of the thermal 
work of the furnace through regulation (stabilization) of the 
temperature of the glass melt in the basic zones of the melt¬ 
ing tank. 
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